Spontaneous Fermi surface deformation in the three-band Hubbard model: A variational Monte 

Carlo study 
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We perform a variational Monte Carlo study on spontaneous c?-wave form Fermi surface deformation ((fFSD) 
within the three-band Hubbard model. It is found that the variational energy of a projected Fermi sea is lowered 
by introducing an anisotropy between the hopping integrals along the x and y directions. Our results show that 
the c/FSD state has the strongest tendency at half-filling in the absence of magnetism, and disappears as the hole 
concentration increases to «/, » 1.15. This is qualitatively in agreement with the mean field analysis and the 
exact diagonalization calculation for the one-band models, and provides a qualitative explanation to the "intra- 
unit-cell" electronic nematicity revealed by the scanning tunneling microscopy. An analysis of the dependence 
of ^/FSD on the parameters of the three-band model indicates that the copper on-site Coulomb interaction, the 
nearest-neighbor copper-oxygen repulsion, and the charge -transfer energy have a remarkable positive effect on 
dFSD. 

PACS number(s):74.20.Mn, 71.10.Li, 71.10.Fd 



I. INTRODUCTION 



Hight temperature superconductivity is realized by doping 
charge caiiiers into an antiferromagnetic insulator In the un- 
derdoped region which lies between hall-filling and optimal 
doping, competing orders occur due to the interplay between 
the localized tendency arising from the strong Coulomb repul- 
sion and the delocalization due to the mobilities of carriers. 
An interesting order tendency is the nematic order,^ which 
breaks the four-fold rotational symmetry of the underlying 
crystal but retains the translational symmetry of the system. 
Recent neutron scattering^ and Nernst effect^ measurements 
of YBaaCusOy provided strong evidences for the existence 
of nematic state in the underdoped cuprates. Moreover, a 
spectroscopic-imaging scanning tunneling microscope (STM) 
investigation on Bi2Sr2CaCu208+a'* suggested an intra-unit- 
cell (lUC) electronic nematic state which breaks the 90° - 
rotational symmetry within every Cu02 unit cell. 

Theoretical studies on the two-dimensional f-/^— and 
Hubbard** models have reported a spontaneous c/FSD in- 
stability in the underdoped region: the Fermi surface (FS) ex- 
pands along the direction and shrinks along the ky direction 
or vice versa. This order is considered to be generated by 
a forward scattering of electrons close to the FS near (0,7r) 
and (tt, 0), which breaks the oiientational symmetry but keeps 
the translational symmetry of the system. From the symme- 
try point of view such a state naturally leads to an electronic 
nematic order However, theoretical studies on the nematic 
order have so far mainly focused on the one-band models or 
in the extreme limit of infinite interactions A recent mean 
field analysis of a more realistic model, i.e. the three-band 
Hubbard modelii, showed that the Coulomb interaction be- 



the delicate mechanism and to depict a clear picture of the 
nematic order in the cuprates. 

To gain a further insight into the nematic order in the 
cuprates, we caiTy out a systematic variational Monte Carlo 
(VMC) study on t/FSD within the three-band Hubbard model. 
Interestingly, we find that the electronic nematicity has al- 
ready shown up at Vpp = 0, and becomes stronger with Vpp 
increasing. We also find that while both the copper on-site 
Coulomb interaction Ud and the charge-transfer energy 
have a remarkable positive effect on t/FSD, the NNN oxygen- 
oxygen hopping tpp has a negative effect. Our results are 
qualitatively in agreement with the slave-boson mean field 
(SBMF) calculations^ over a wide doping range, showing 
that (iFSD has the strongest tendency around the van Hove 
Ijjjjjjg8j9ji2,i6ji7 (namely half-filing in our calculations) in the 
absence of magnetism. Beyond that, we show that dFSD in 
the three-band model leads to an imbalance in the hole den- 
sities of the neighboring oxygen sites, providing a consistent 
theoretical explanation on the phenomena observed by STM 
in Bi2Sr2CaCu208+A 

The paper is organized as follows: In Sec. II we define the 
three-band Hubbard model and outline the VMC scheme. Re- 
sults from this numerical solution and comparisons with pre- 
vious studies are presented in Sec. III. A conclusion is given in 
Sec.IV. 



tween next nearest-neighbor (NNN) oxygen atoms Vpp plays 
a crucial role in the formation of the nematic order , which 



demonstrates that the interaction within the unit cell on the 
Cu02 plan is needed to be taken into account in order to show 



II. MODEL HAMILTONIAN AND METHODS 



By considering the energetics and hybridizations for copper 
3t/j.2_,,2 orbital and oxygen Ip^ (2py) orbital, the kinetic part 
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of the three-band model reads, 

i,(T V 

-hp 2 Yj {ply/2,aP>+y'/2.<r + h-^), (1) 
ucr {vy) 

where d' and p.' create a hole with spin cr at the /th 
copper site and the / + v/2th oxygen site, respectively, with 
V = flx(ay) being the unit vectors along the x and y direc- 
tions, n'l^ and n'!^y/2o- corresponding number opera- 
tors. tp4 and tpp denote the copper-oxygen and oxygen-oxygen 
hopping integrals, (v, v') limits the sum over NNN oxygen- 
oxygen lattice sites. We set = on the copper site and 
introduce the charge-transfer energy A^, - Sp - Sd to control 
the relative Cu/O hole densities. 

The interaction part of the model is given by, 

H' = UdY ntn' + ^ y < „ , 

/ i,v 
i,v o-,a' 

+ ypp y y n'' n1 „ (2) 

/ i / j (+v/2,(T /+y/2,cr" ^ ' 

here f/^and f/p are the on-site Coulomb interactions on the 
copper and oxygen sites, respectively. Vpd and Vpp are the 
Coulomb repulsions between nearest-neighbor (NN) copper- 
oxygen sites and NNN oxygen-oxygen sites, respectively. In 
the following, the parameters within CuOa plane are measured 
in the unit of tpd- 

We use the following variational wave function: 

\k\<kF,(T 

where Pq - PdPpPpdPpp is the projection operator, which 
contains four parts: 1. Pd - FI (l ~ (1 ~ ^d) nf|nf|) mod- 

ifies the double occupancy on the Cu site; 2. Pp - 

n(l - (l -^p)nf+^/2.t«f+v/2,i) modifies the double occu- 

d P 

pancy on the O site; 3. Ppd = Yl gpd"'"'^''^- adjusts the 
charge occupancies on the NN copper-oxygen sites; 4. Ppp = 
n gpp"'*''^"'*'''^ adjusts the charge occupancies on the NNN 

i,<v,v'> 

oxygen-oxygen sites. Here gd, gp, gpd, and gpp stand for four 
variational parameters in the range from to 1, which will be 
optimized during the Monte Carlo calculations. 

In Eq. (O, the single-particle operator o-^^ is given by a lin- 
ear combination of dj^^, p[i^^, and p'^j^^, which are the Fourier 
transformations of d- and p- hole operators: 



where is the total number of unit cells. The coefficients in 
ff^^ are obtained through diagonalizing the following matrix, 

Aeff ^ -4fpp COS |- COS Y -2fprf,.v COS 1^ ' 

-4tpp COS |- COS -J Aeff -2tpd,y COS -j- , 

, -2fp^_.i. COS |- -2tpd,y COS ^ 

(6) 

here A^// is a variational parameter acting as an effective 
charge-transfer energy. Introducing an asymmetry parameter 
5var between the hopping integrals along x and y direction, we 
have 

tpd,x - tpd - <5var, (7) 
tpd,y - tpd + 5var- (8) 

When 6yar is finite, the shape of FS is deformed. We evalu- 
ated the expectation value of total energy for the variational 
wave-function by employing Monte Carlo method The 
four projection parameters gd, gp, gpd, and gpp, as well as the 
effective charge-transfer energy A^ff and the asymmetry pa- 
rameter dyai; are optimized to obtain the minimal energy. Dur- 
ing the optimization, a quasi-Newton method combined with 
the fixed sampling method"^ -^ is used. The calculations have 
been done for square lattices with periodic and antiperiodic 
boundary conditions along the x and y direction, respectively. 
These boundary conditions are used to avoid the degenerate 
states at the Fermi surface. Typically, 10^ Monte Carlo steps 
are performed for each set of variational parameters. The re- 
sulting statistical errors are given in the relevant figures by 
error bars. 



III. NUMEmCAL RESULTS AND DISCUSSIONS 

In the present paper, we take the typical values of the copper 
oxide parameters. Setting the copper-oxygen hopping = 1, 
we have Ud = 8.0, Up = 3.0, Vpd = 1.0, and A„ = 3.0, 
according to the constrained density-functional calculations^^. 
tpp is expected to have a negative effect on c/FSD, so we will 
firstly focus on the case of f^p^O.O and the order dependence 
on tpp will be examined alone later on. The results presented 
below are obtained for these typical parameters, except where 
explicitly noted otherwise. 

Firstly, in Fig[T] (a), we present the doping depen- 
dence of the condensation energy per unit cell {Econd - 

((5vfi') - £(0)j IN, with (Syfi' being the optimized asymme- 
try parameter) on the square lattices of the sizes of 16x16, 
20x20, and 22x22. One can observe that for all the three lat- 
tices, the absolute value of Econd decreases continuously from 
half-filling and vanishes at «/, ~ 1.15. This behavior is con- 
sistent with the slave-boson mean-field (SBMF) analysis for 
the t- J model, indicating that the growing forward scattering 
with decreasing hole density has a positive contribution to the 
formation of the c/FSD state. Fig[T](b) shows the dependence 
of (JX' on the hole density. There is an obvious tendency that 
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FIG. 1 : (Color online) Condensation energy Ecomi (a) and optimized 
value of (5viu (b) as a function of hole density on the 16x16, 20x20, 
and 22x22 lattices. 



(JX' decreases to zero at ~ 1.15 in the same way as the 
condensation energy, suggesting that c/FSD only occurs in the 
underdoped region. 
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FIG. 2: (Color online) Fermi surfaces in the normal state ((5vai =0.0) 
(a) and the dFSD state (5v„r=0.062) (b) at the hole density «,,=1.05. 
Filled (empty) sites indicate the occupied (unoccupied) k points. 

We present the isotropic (5vai=0.0) and distorted 
(5y^;.=0.062) FS's at nh-l.05 for a 20x20 square lattice 
in Fig|2] It is shown that the four fold rotational symmetry 
of FS is reduced to two fold in the distorted ground state. 
The FS expands along the direction and shrinks along 
the direction as the optimized 5var takes a positive value. 
We learn from Figi2]that the finite size effect brings some 
difficulties in determining the optimal value of 5vai in our 
calculations: when the FS is distorted by a nonzero 5vai, the 
states near (0,7r) or (7r,0) undergo a discontinuous change 
from unoccupied to occupied or vice versa. This kind of 
discontinuity introduces an uncertainty of the optimal value 
of (5var, which is estimated about 0.01 ~0. 02 for our studied 
systems. 

From FiglT] it can be seen that the c/FSD state is stablest 
at half-filling. However, strong (n, n) scattering around half- 
filling may favor an antiferromagnetic (AFM) state, instead of 
the (iFSD state. To clarify this issue, we carried out VMC cal- 



culations for the projected AFM state at half-filling, and the 
corresponding condensation energy per unit cell is -0.2185 
and -0.2266 on the 16x16 and 20x20 lattices, respectively. 
Two order lower of the condensation energy in the AFM state 
compared to the t/FSD state demonstrates that magnetic or- 
dering is actually the strongest instability close to half-filling. 
With the increase of doping concentration, the AFM order 
quickly disappears in doped cuprates, and the ground state 
could be well described by a wave function without long- 
range magnetic ordering, such as the one in Eq. (|3). 

Another long-range order that may compete with t/FSD is 
the (/-wave superconductivity, which is induced by AFM spin 
fluctuations. According to the VMC study of the one-band t-J 
moder, the c/FSD instability is overwhelmed by the c/-wave 
superconductivity. However, in the case of the three-band 
Hubbard model, a VMC study^^ showed that the supercon- 
ducting condensation energy is about 10"^, which is compa- 
rable to the energy gain in the c/FSD state by our calculations. 
This indicates that the c/FSD instability is as strong as the su- 
perconducting instability in the three-band Hubbard model. 
A further study is needed to clarify one remaining question, 
whether both c/FSD and c/-wave superconductivity can coex- 
ist in the multi-band model. 

To understand the physical origin for the formation of 
c/FSD, we present different energy changes AEa - (^^(^vai) - 
Ea(0))/N as a function of <5var, with a representing different 
components of the Hamiltonian in Fig. [3] We observe that 
while AEJcin and hEJJp take positive values as the system 
changes into an anisotropic state, the other three components 
AE_Ud, AEJVpd, and AE_Ae, take negative values in the stud- 
ied parameter region, which lower the total energy of the sys- 
tem, and lead to a minimum at a certain value of dym - This re- 
veals that the c/FSD state is the result of delicate competition 
between the localized tendency introduced by Ud, Vpd, and 
the delocalization due to the kinetic term of carriers. One can 
notice that the dependence of different energy changes on dyai 
shows similar behavior both at half-filling (n;, = 1) and at a 
finite hole doping density (n/, = 1.05). This is rather different 
from the VMC study of the t-J model at half-filling^, where 
the kinetic energy is exact zero and the J induced superex- 
change energy is found to increase with <5vai increasing, re- 
sulting in the vanishing of condensation energy at half-filling. 

To explore the parameter dependence of the c/FSD state, we 
present the change of total energy AE - [E (5var) - E(0)] /N 
as a function of (5vai for different sets of parameters at a fixed 
doping («/, = 1.05) in Figs. |4l|6] From FigH] one sees that 
the on-site Coulomb interaction on the copper sites, Ud, has a 
remarkable positive effect on c/FSD, manifested by a decrease 
of AE with increasing Ud and a shift of minimum to larger 
dyai - This can be expected since Ud makes a dominant con- 
tribution to the forward scattering, which is considered to be 
the origin of c/FSD. As seen in Fig.|5] a decrease of AE with 
increasing Vpd or A^, indicates that both NN copper-oxygen 
repulsion and charge-transfer energy have a positive effect on 
the formation of c/FSD. It is naturally to expect that when Vpd 
or Aci increases, more holes transfer from oxygen sites to cop- 
per sites, resulting in a stronger forward scattering of carriers 
hyUd. 
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FIG. 3: (Color online) Different energy changes as a function of Sy^y 
on the 20x20 lattice, (a) ;j;, = l, (b) n;, = 1.05. 




FIG. 4: (Color online) Total energy change as a function of (5var for 
different values of C/^. Lattice = 20x20 and n;,=1.05. 



Once tpp turns on, it may enhance the delocalization of 
holes on the oxygen sites, and have a strong negative effect on 
t/FSD. To examine if t/FSD exists in the actual cuprates mate- 
rials, we compare the change of total energy AE for tpp - 0.0 
and tpp = 0.4 in Fig. |6] It is readily seen that although the 
condensation energy at tpp-QA is much smaller than the case 
of tpp-0.0, the c/FSD state still survives. 

The NNN oxygen-oxygen Coulomb repulsion Vpp is ex- 
pected to have a positive effect on t/FSD, since it favors an 
inhomogeneity in the hole densities of the NNN oxygen sites. 
In the previous work by Fischer et al.'^, a mean field analy- 
sis based on the three-band model showed that the minimum 
value of Vpp that needed for a system to enter a nematic phase 
is about 1.2, which is much higher than the realistic value. 
Here, in Fig|6] we show that the electronic nematicity has 
already shown up at V^^^O, and becomes stronger if we in- 
troduce Vpp-Q.5. Our results demonstrate that a spontaneous 
t/FSD can occur in actual cuprates materials. 




FIG. 5: (Color online) Dependence of total energy change on 
for different values of V,,j (a) and A„ (b) on the 20x20 lattice and 

)l;,= 1.05. 
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FIG. 6: (Color online) Total energy change as a function of Sy. 

20x20 and «i = 1.07. 



for 



different values of tpp and Vpp. Lattice 



As mentioned before, an STM investigation on 
Bi2Sr2CaCu208+a^ showed that there exists a significant 
density difference in the two oxygen sites within each unit 
cell. In Fig|2l we present the hole densities on the copper 
and oxygen sites for the 20x20 lattice, at n/,=1.05. The 
inhomogeneous distribution of hole densities in the p ^ and 
oxygen orbitals indicates that liFSD in the three-band 
model naturally results in an lUC nematicity, which is in 
good agreement with the results of the STM investigation. In 
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the investigation of high-Tc cuprates, nematic order is usually 
considered to be the result of melted stripesM, namely, 
arising from proliferation of dislocations, the topological 
defects in the striped state^*^. However, although a striped 
state may have a close relation to the nematic order, the t/FSD 
instability can also lead directly to an electronic nematic 
state, as shown in Fig|7] 



X 

FIG. 7: (Color online) Distribution of hole densities on the copper 
and oxygen sites. The results are obtained for the 20x20 lattice at 
«;,=1.05. The length of blue squares represents the hole density at 
O, and Ov sites, red circles indicate the locations of Cu sites. 



parameter region of cuprates, which is weakened with increas- 
ing the hole doping density and vanishes at «/, ~ 1.15. The re- 
sulting imbalance of hole densities on the NNN oxygen sites 
confirms that t/FSD is a possible origin for the lUC nematic 
phenomena, which was firstly proposed by Lawler et al^. Our 
results show a qualitative agreement with the experimental re- 
sults^^ and are consistent with the SBMF calculations for the 
t-J model''. An analysis of the order dependence of liFSD on 
various parameters within the three-band model showed that 
Ud, Vpd, and A^., have a remarkable positive effect on liFSD. 
They contribute to the formation of t/FSD in such a way that 
the forward scattering is enhanced either by increasing Ud or 
through transferring holes from oxygen to copper sites as Vpd 
(Act) becomes larger We found that c/FSD survives when a 
finite oxygen-oxygen hopping integrals tpp is introduced, and 
if a finite Vpp interaction exists, the c/FSD state will be stabler 

Recently, the experimental evidence for electronic 
anisotropy in iron-based superconductors has been accu- 
mulated^'^-''*. The situation in iron-based superconductors 
is much more complicate since the lattice, spin, and orbital 
degrees of freedom all manifest themselves in the nematic 
stated. We consider that our VMC study for the cuprates 
can be generalized to iron-based superconductors and may be 
helpful for clarifying the nematicity in these systems. 
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